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The ground state of a spin i nearest neighbor quantum Heisenberg antiferromagnet on the py- 
rochlore lattice is investigated using a large A'^ SU (N) fermionic mean field theory. We find several 
mean field states, of which the state of lowest energy upon Gutzwiller projection, is a parity and time 
reversal breaking chiral phase with a unit monopole flux exiting each tetrahedron. This "monopole 
flux" state has a Fermi surface consisting of 4 lines intersecting at a point. At mean field the low- 
energy excitations about the Fermi surface are gapless spinous. An analysis using the projective 
symmetry group of this state suggests that the state is stable to small fiuctuations which neither 
induce a gap, nor alter the unusual Fermi surface. 



I. INTRODUCTION 



This paper lies at the intersection of two streams of re- 
search in contemporary quantum magnetism — the study 
of spin liquids and the study of geometrically frustrated 
magnetism. Specifically, we are interested in 5 = 1/2 
Heisenberg models on the pyrochlore lattice and were 
motivated by asking whether they support a zero temper- 
ature phase that breaks no symmetries of the problem — 
a fully symmetric quantum spin liquid. 

The study of quantum spin liquids — being defined 
broadly as states of spin systems that do not exhibit long 
range Neel order down to zero temperature — iss currently 
in the midst of a significant revival. The subject itself is 
decades old with its contemporary study tracing its ori- 
gins to Anderson's introduction of the resonating valence 
bond (RVB) stated and then to his suggestion^, upon the 
discovery of the cuprate superconductors, that their be- 
havior was traceable to a parent spin liquid state. But its 
current vogue has much to do with recent progress in con- 
structing actual models that realize spin liquid behavior^ 
and the recognition that a large class of spin liquids exem- 
plify ordering beyond the broken symmetry paradigm — 
they give rise to low energy gauge fields but not order 
parameters. That such "topological phases" ^ also un- 
derlie a fascinating approach to quantum computation^ 
only multiplies their interest. 



^ For an introduction to this area see the recent Les Houches lec- 
tures by Misguich^ and the older review article Ref. 0. 

^ Here we use the term "topological phases" in the looser sense of 
any phase with emergent gauge fields. Strictly speaking the term 
should be reserved for cases where the low energy gauge theory 
is a purely topological gauge theory. 



The study of geometrically frustrated magnets^ has in- 
tertwined roots. Indeed, Anderson's 1972 paper identi- 
fied a small value of the spin and geometric frustration 
as two sources of quantum fluctuations that could fa- 
vor a spin liquid. In recent years there has been steady 
progress in both understanding the behavior of many ge- 
ometrically frustrated magnets but, more importantly, 
in synthesizing an increasing number of compounds that 
realize challenging idealizations to increasing accuracj*2i^ 
leading to a resurgence of interest in these systems as 
well. 

The pyrochlore lattice is a natural object of study in 
this context. It is highly frustrated and frequently re- 
alized as a sublattice of the spinels or the pyrochlores. 
Potentially, it could host a spin liquid in c? = 3 for small 
values of the spin. Much work has gone into studying its 
magnetic properties in various contexts. Most notably, 
it is known to lack long range order with nearest neigh- 
bor interacting classical spins^ but instead to exhibit an 
emergent gauge field and dipolar correlations as T 0. 
[Interestingly, this physics is realized in the Ising "spin 
ice" systems (Dy and Ho titanate)^ although with an 
additional fundamental dipolar interaction that leads to 
further elegant physics involving magnetic monopolesii.] 
Attempts to work about the classical limit, in the spin 
wave (l/S") expansion have lead to some insight into 
the quantum "order by disorder" selection mechanism 
in this limit. While the fate of the 1/S expansion is not 
scttlcdii^, there is however little reason to think that it 
can be informative when it comes to small values of spin, 
especially the 5=1/2 case'^ that is our concern in this 



^ This is as good a place as any to note that there is not, to date, 
a good experimental S = 1/2 antiferromagnet on the pyrochlore 
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paper. This is so partly because the selection mechanism 
at large S is weak and leads to somewhat ornate states 
but also for the well-understood reason that it misses out 
on tunneling processes that are sensitive to the Berry 
phases entering the exact path integra l^^i^^ . 

Consequently, various authors have attempted to di- 
rectly tackle the 5=1/2 problem. Harris. Berlinsky and 
Bruder— initiated a cluster treatment in which the py- 
rochlore lattice is first decoupled into, say, its up tetrahe- 
dra and then perturbatively reconnected. Subsequently 
Tsunetsugrti^ worked out a more complete treatment 
along the same lines and found a dimerized state with 
a four sublattice structure. The criticism that this work 
predicts symmetry breaking that is put in at the first step 
has attracted a potential rebuttal in the work of Berg 
et alii with the "Contractor Renormalization" or CORE 
technique. An alternative perspective on this physics was 
provided \n-^ where it was shown that an SU{N) defor- 
mation produces a quantum dimer model whose physics 
is very reminiscent of the HBB scenario. Unfortunately, 
the N = 2 limit is manifestly problematic so it has not 
been possible to declare victory in this work. Yet an- 
other attack on the probleroi^ used an alternative large N 
theory — equivalent to Schwinger boson mean field theory 
— and found a delicate energetics at small values of spin 
(or boson density) which nevertheless strongly indicated 
that the spin 1/2 problem must break some symmetry.^ 

With this set of predictions of symmetry break- 
ing as background, in the present work we bring an- 
other approximate large N technique — that of "slave 
fermions"— — to bear on the pyrochlore problem with 
a view to examining whether it produces a symmetric, 
spin liquid, alternative. To this end we enumerate vari- 
ous translationally invariant mean field solutions of which 
the lowest energy non-dimerized solution is one we call 
a "monopole flux" state; upon Gutzwiller projection it 
also improves upon the fully dimerized states. While 
this state does not break lattice symmetries in the man- 
ner of the HBB scenario, it is not a spin liquid in the 
sense of breaking no symmetries at all. Instead it is a 
chiral spin liquid^^i^ and breaks parity (P) and time re- 
versal (T) symmetries. It also exhibits spinous in its 
mean field spectrum. We describe the unusual mean 
field spectrum — which yields a Fermi surface consisting 
of four lines intersecting at a point — and its low energy 
limit in some detail. This state was first reported in 
Ref. [l^. Subsequently it sparked a larger investigation 
by R. Shankar and two of us^^ on flux Hamiltonians on 
root lattices of Lie groups with minuscule decorations 
and these results were announced there previously. The 
stability of the mean field structure to fluctuations is the 



lattice. 

Ref. [l^ shows that at asymptotically small boson densities the 
system must break some symmetry. The minor caveat is that 
this does not rule out a different solution intervening right near 
S = 1/2. 



next question of interest. We make progress in that di- 
rection by enumerating the projective symmetry group 
(PSG)^ of the state and showing that it forbids any 
terms that would destabilize the mean field Fermi sur- 
face. This still leaves the fate of the gauge fluctuations 
open as a matter of dynamics and we expect to discuss 
this elsewhere^l. Finally we note that as we were finaliz- 
ing this paper there appeared an independent evaluation 
of energies for Gutzwiller projected wavefunctions on the 
pyrochlore^^ which agrees with our results on that score. 

And now to the organization of the paper. We be- 
gin with a brief overview of the large-N/mean field slave 
fermion treatment of the Heisenberg model in Section [Til 
In Section [mi we apply this technique to generate several 
mean field ansatze on the pyrochlore lattice. We iden- 
tify the lowest energy state, or monopole flux state, and 
discuss its interesting properties. Section IIVI reviews in 
general terms how the PSG protects a mean fleld state 
against developing symmetry-breaking terms. The PSG 
derived arguments for the stability of the monopole flux 
state are given in Section |Vl where we derive the gen- 
eral form of the symmetry permitted perturbations to 
the Hamiltonian. We conclude in Section IVIl Details 
of the PSG for the monopole flux state can be found in 
Appendix |B] while Appendix [A] explains the numerical 
technique used to carry out Gutzwiller projection. 



II. THE LARGE-N HEISENBERG MODEL: 
SPINONS AND GAUGE FIELDS 

In this section we briefly review the large N fermionic 
approach to the 5=1/2 SU{2) Heisenberg model which 
began as a mean field theory introduced by Baskaran, 
Zou and Anderson^S and was shortly thereafter system- 
atized via a generalization to SU{N) by Affleck and 
Marston^i. 

In this approach, we first replace the bosonic spin op- 
erators of the Heisenberg Hamiltonian 

= J ^ S, . Sj (1) 

with bilinears in fermionic "spinon" operators: 

Sz = i ^ 4^Crai3Ci,3 ■ (2) 
a, (3 

The resulting Hamiltonian conserves the number of 
fermions at each site and the starting spin Hamiltonian is 
recovered if we limit ourselves to physical states with ex- 
actly 1 particle per site. Up to a constant in the subspace 
of physical states, it can be re- written in the suggestive 
form, 

^ = X] 4aCiac]^C,^( (3) 

<ij> a p 
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A mean field theory arises upon performing the Hubbard- 
Stratonovich decoupling 



a <i]> 



<ij> 



and locally minimizing the classical field Xij to obtain 
self-consistency. 

In order to understand the nature of fluctuations about 
such mean field solutions it is conceptually convenient 
to consider the path integral defined by the equivalent 
Lagrangian: 

L = ^4„(0t +/^)c,.„ +^0j(4„c,,„ - 1) 



E 

<ij> 



J 



(5) 



where is a Lagrange multiplier field enforcing the single 
occupancy constraint J2a Aa^ia = 1- 

The above Lagrangian ([5]) is invariant under the local 
gauge transformations 
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(t>~* 4> + de/dt 



(6) 



which arise from the local constraints in the fermionic 
formulation. It follows that we have reformulated the 
Heisenberg model as a problem of fermions that live on 
the sites of the original lattice coupled to a C/(l) gauge 
field and an amplitude field (the phase and amplitude of 
Xij) that both live on the links of the lattice. In other 
words, we may write Xij = Pij^'''^^^, where aij + 

Oi — 6j under the gauge transformation ([6|. The mean 
field theory consists of searching for a saddle point with 
frozen link fields. 

As the Lagrangian ([5]) docs not directly constrain the 
phase of the Xiji it describes a strongly coupled gauge 
theory where the assumption of a weakly fluctuating 
gauge field invoked in the mean field theory is, prima 
facie, suspect. To circumvent this barrier, Affleck and 
Marston^i proposed a large N framework which intro- 
duces a weak-coupling limit for the model ^ by extend- 
ing the SU(2) spin symmetry group of the Heisenberg 
model to SU(N) with N even. The result is a theory 
of many spin flavors whose coupling strength scales as 
J — > J /N . In the limit that iV — *■ oo, the correspond- 
ing mean field theory is exact; for sufficiently large but 
finite N one hopes that a perturbative expansion gives 
accurate results. The validity of the qualitative features 
deduced at large N in the starting SU{2) problem is, 
of course, hard to establish by such considerations and 
requires direct numerical or experimental confirmation. 

To effect the large N generalization, we replace the 2 
spinon operators C| and C| with N spinon operators Cq. 



The single occupancy constraint is now modified to 

N 

2 



(7) 



and the large-N Hamiltonian has the form 

Q,/3 <ij> 



= - E E (4cjax»j- + h.c.) + 1' 



a <ij> 



J 



<ij> 



In the infinite N limit, the action is constrained to its 
saddle point and the mean field solution becomes exact. 
Further, to lowest order in the allowed fluctuations 
involve moving single spinous, so that as iV — *■ oo we 
need only impose the constraint ([7]) on average. 

Away from N ~ oo the link fields, especially the gauge 
field, can fluctuate again although now with a control- 
lably small coupling. While the fate of the coupled 
fcrmion-gauge system still needs investigation, the pres- 
ence of a small parameter is a great aid in the analysis, 
as in the recent work on algebraic spin liquids2^. 

Finally, we note that the starting SU{2) problem is 
special, in that it is naturally formulated as an SU{2) 
gauge theory^i^. This can have the consequence that 
the N = 2 descendant of the large N state, if stable, may 
exhibit a weakly fluctuating SU(2) gauge field instead of 
the U{1) field that arises in the above description. We 
will comment on this in the context of this paper at the 
end. 



III. MEAN-FIELD ANALYSIS 

A. Saddle Points of the Nearest Neighbor 
Heisenberg Model 

We begin by enumerating mean field (MF) states which 
preserve translation invariance on the pyrochlore. A 
mean field solution consists of a choice of link fields which 
minimizes the mean field energy functional for the La- 
grangian ([5]) 



Eix) = N 



Ejl^yl 



E(^(k)-/^) 



(9) 



where e(k) is the energy of a spinon of momentum k in 
the fixed background Xiji ^^'^ the chemical potential fi 
is chosen so that the constraint of 1 particle per site is 
satisfied on average. 

As discussed in the introduction, previous work on the 
Heisenberg model on the pyrochlore lattice has led to 
ground states with broken symmetries. In this work we 
are particularly interested in constructing a natural state 
on the pyrochlore that breaks as few symmetries as pos- 
sible. To this end, we begin our search with especially 
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symmetric ansatze for which pij = p is independent of 
i and j, and the flux <i>A = X^a "^y through each face 

of the tetrahedron is the same. The net flux J^Ui *A 
through each tetrahedron must be an integer multiple of 
2tt, since each edge borders two faces such that its net 
contribution to the flux is (mod 27r). This gives the 
following 3 candidate spin liquid states: 



1. Uniform: <f> 



A 







2. TT Flux: $ 



A 



3. Monopole: $a = Tr/2. Every triangular face of the 
tetrahedron has a 7r/2 outwards flux - equivalent 
to a monopole of strength 2tt placed at the center 
of each tetrahedron. 

At infinite N a dimerized state is always the global min- 
imum of thus we also consider 

4. Dimerized: Xij = on a, set of bonds that con- 
stitute any dimer covering of the lattice but zero 
otherwise. 

The states (1-3) above are analogues of the uniform, tt 
flux, and chiral states studied previously on the square 
Iattice2ii2^. Of the above states, (1) and (2) break no 
symmetries of the problem; the third preserves lattice 
symmetries but breaks P and T. 

The states (1) and (2) arc in fact particle-hole conju- 
gates: a particle-hole transformation maps clcj + CjCi — > 

—CjCi — clcj, changing the sign of x on each bond and 
adding tt flux to each triangular plaquette. At = 2 this 
can be effectuated by an SU (2) gauge transformation, so 
that the states (1) and (2) describe the same state after 
Gutzwiller projection. 

The mean field energies of these states are listed in the 
first column of Table I. Consistent with Rokhsar's gen- 
eral considcrations^^i^ the fully dimerized state is lowest 
energy and the monopole flux state has the lowest en- 
ergy of the non-dimerized states. The mean field states 
with N set equal to 2 do not satisfy the single occupancy 
constraint. While, in principle, perturbation theory in 
1 /N can greatly improve the wavefunction in this regard 
this is a complex business (to which we return in Sec- 
tions IV and V) ill-suited to actual energetics. Instead, 
the somewhat ad hoc procedure of (Gutzwiller) project- 
ing the mean field wave function onto the Hilbert space 
of singly occupied sites is typically employed to improve 
matters. This leads to resonances and long range correla- 
tions that can substantially lower the mean field ground 
state energy, particularly for spin-liquid type states. 

Expectation values in the Gutzwiller projection of a 
state can be carried out using a Monte Carlo approach, as 
described in Ref. lU. A brief description of the numerical 
method specialized to our problem is given in Appendix 
1X1 The second column in Table U shows the numeri- 
cally evaluated energies of the 4 mean field states with 
Gutzwiller projection. We see that the monopole flux 





Emf (unprojected) 


Emf (projected) 


Uniform 


-0.3333J 


-.3752 ±0.0004 


TT Flux 


-0.3333J 


-.3752 ±0.0004 


Monopole 


-0.3550J 


-.4473 ±0.0009 


Dimer 


-.375J 


-.375J 


(7r,7r) 


-0.3333J 


-0.3751 ± 0.0008 


{n/2,7T) 


-0.3491J 


-0.4353 ±0.002 



TABLE I: Mean-field energies for projected and unprojected 
ground states of the mean field ansatze considered. The 
quoted mean field energies are the energy of (j4| plus the 
omitted constant per site required to make a correspon- 
dence with ([T}. The states (tt, tt) and (7r/2,7r) are variants 
of the uniform and monopole state, respectively, with fiux n 
per hexagonal plaquette. The projected wave functions were 
evaluated on a lattice of 5 x 5 x 5 unit cells, or 500 sites 
for configurations with a 4 site unit cell, and 1000 sites for 
configurations with an 8 site unit cell. 



state now emerges as the lowest energy state of our quar- 
tet. Encouraged by this, and also because the state has 
various elegant properties, we will focus in the remain- 
der of this work on the properties of the monopole flux 
state. Note however, that we have failed to preserve all 
symmetries of the Hamiltonian even in this approach — 
we are forced to break T and P and thus end up with 
a chiral spin liquid. We give a fuller description of the 
symmetries of the state below. 

Finally, we note that larger unit cells can be consistent 
with translationally invariant states.^ Such states have 
an integral multiple of 7r/2 flux through each triangular 
plaquette, but also non-trivial flux through the hexagonal 
plaquettcs in the kagomc planes, as for the mean field 
states on the kagomc studied in Ref. [s^. By the same 
arguments as employed for a single tetrahedron, we find 
that the flux through the hexagons must have values or 
TT (mod 27r) to preserve the translational symmetry of the 
lattice. [A flux of tt/2 per hexagonal plaquette necessarily 
breaks lattice translations]. However, as noted in Tabled 
wc find that these states also have higher energies than 
the monopole flux state both at mean fleld and upon 
Gutzwiller projection. 



B. The Monopole Flux State 

The monopole fiux state exhibits a flux of 7r/2 per tri- 
angular face. To write down the mean field Hamiltonian 
explicitly we must pick a gauge. We choose Xij = Poe'"'^ , 
with aij = ±7r/2. The phase of ±i that a spinon picks up 
when hopping from site i to j can be represented as an 
arrow on the corresponding edge, which points from i to 



We are grateful to Michael Hermele for emphasizing this point. 
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FIG. 1: Link field orientations in the monopole state with 
X — ±ipo- Hopping along the direction of an arrow induces a 
phase of n/2; hopping against the arrows, a phase of —-k 12. 
The flux on each triangular face is 7r/2 outwards. With this 
flux assignment the monopole flux state breaks T and P, but 
is invariant under lattice translations and rotations. 



j [j to i) if the resulting phase is +(— )i. The orientation 
of the link fields, shown in Figure [1] gives an outward 
flux of 7r/2 per plaquette. 

The necessity of picking a gauge for the mean field so- 
lution causes, as usual, various symmetries to be imple- 
mented projectively. For example, the assignment shown 
in Figure [T] is not invariant under lattice rotations. How- 
ever, the background link fields after rotation can be 
gauge transformed to the original state, as expected from 
the manifestly rotation invariant assignment of fluxes. 
We discuss these and other symmetries in more detail in 
Sections II VI and fVl here we merely note that P and T are 
the only symmetries broken by the monopole flux state. 

The Hamiltonian for spinous in the gauge choice shown 
in Fig. [l]is 



2Np 





sin( '°"1^^^ ) sin 



sin(^) sin(^) sin(^) 
sin(^) sin(^) 
sin(^^) 



ka • 



(10) 



sin(4^) sin(^^) sin(^) 







where is a 4-component vector, with = c^^. Here 
the index i labels the 4 sites in the tetrahedral unit cell. 



Low Energy Expansions of the Spinon 
Dispersion 



Figure 2(b) shows a plot of the energy eigenvalues of 
(I10[) along the high-symmetry lines of the Brillouin zone. 
At half filling, the Fermi 'surface' consists of the lines 
fc(±l,±l,±l) which join the point (0,0,0) to the center 
of the hexagonal faces of the Brillouin zone of the cubic 
FCC lattice, line (L - F) in Fig. |2(b)[ Each Fermi hue 
has a pair of zero energy eigenstates. 



The low-energy structure of the monopole flux state 
can be divided into two regions: the set of 4 Fermi 
lines sufficiently far from the origin, and i?o, the area 
near the origin. 

In only 2 of the four bands lie near the Fermi 
surface, and the low-energy theory is effectively two- 
dimensional. Linearizing the Hamiltonian about one of 
the Fermi lines gives: 



Figures 3(a) shows a surface of constant energy i? w 
near the Fermi surface. At = 0, the 4 bands inter- 
sect only at the origin and the constant energy surface 
is given by the 4 lines described above. Surfaces of con- 
stant energy E k Q^E ^ consist of 4 cylinders enclosing 
the (1, 1, 1) directions, which are the surfaces of constant 
energy for particle-like {E > 0) or hole-like {E < 0) ex- 
citations about the Fermi line. About the origin all 4 
bands have energy linear in fc, and another, diamond- 
shaped constant-energy surface appears. These surfaces 
intersect at the band crossings along the x, y, and z axes^. 



This Fermi surface does not display fermion doubling in the 
naive sense; all four bands cross at only one point in the Brillouin 
zone. This does not violate the result of Ref. |36L which assumes 
that levels are degenerate only at a finite set of points. 



H[s,e] = [fc, e, 6*] (e cos 6Iti -he sin 6'r2)*ia[fc,e, 61], 

(11) 

with energies ±e, independent of 6. Here we have used 
the local coordinate system 



/ 2 u 1) u V 

with 9{u,v) = tanr^{v/u) and e(u,w) = 
\/u^ + I (2>/2). Curiously, at mean field the low 
energy spectrum is independent of the position k along 
the Fermi line, depending only on the momentum com- 
ponent in the kagome planes perpendicular to the vector 
Zj. Thus the linearized theory away from the origin 
consists of a continuum of flavors of Dirac fermions 
confined to the kagome planes orthogonal to this line. 
In Rq all 4 bands have energies vanishing linearly as 




FIG. 2: (b) shows the spectrum of the monopole flux state. 
Note the Fermi hue -^{k, k, k). (a) shows the contour in the 
BriUouin zone along which the spectrum is plotted. 



fc ^ 0, and the low-energy Hamiltonian is given by: 



H 



2Np 





4 

, 4 



4 





k 





k ^k 
4 





with energy eigenvalues 



e = ± 



(13) 

(14) 



(»<j) 



This dispersion relation also gives massless spinous; how- 
ever, the theory is no longer one of Dirac fermions. 

In addition to four bands touching at the origin, the 
linearized Hamiltonian (|13|) has 2 zero eigenvalues on 
each Fermi line. Restricting the spinors to the corre- 
sponding low-energy subspace again yields the expansion 




(a) 




(b) 



FIG. 3: Constant energy surfaces of the monopole flux state, 
for E/ J = 0.5. Decreasing E/ J makes the cylinders thinner, 
(a) A view of the cube of side length tt surrounding the origin. 
Eight cylinders surrounding the eight Fermi lines emanate 
from the origin; at the origin a diamond-line shape (the low- 
energy spectrum of the remaining two bands) can also be seen. 
This shape repeats at the cubes corners (ivr, ±7r) . (b) A close- 
up view of the region surrounding the origin. Altering the 
ratio E/ J shrinks the entire structure, but does not change 
its shape. 



(fTTjl . Thus captures the principal features of the low- 
energy behavior not only in the vicinity of the origin, but 
throughout the entire BriUouin zone. 

The linearized Hamiltonian has several interesting fea- 
tures. First, we may express it in terms of three matrices 
as follows: 



H = Uxkx + ayky + a^kz 



(15) 



The a matrices arc reminiscent of Dirac 7 matrices, albeit 
with a tetrahcdrally invariant, rather than rotationally 
invariant, algebraic structure. They do not comprise a 
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Clifford algebra, but obey the anti-commutation relations 



24 



V3\e,,k\Wk 



(16) 



Further, in a 3+1 dimensional Dirac theory there are 



two matrices (70 and 75) which anti-commute with all 7^. 
In this sense our mean field Hamiltonian more resembles 
a 2+1 dimensional Dirac theory: there is a unique matrix 
ao such that {ao,ai} = 0,i = 1..3, given by 



v(0) 



cos(y) cos(a;) 

cos{y) cos(a;) 

- cos{y) cos(z) — cos(a;) 003(2;) 

' cos(a;) cos(z) cos(j/) 003(2:) 



cos(j/) 003(2) 
cos(a;) 003(2) 


— C03(y) C03(x) 



cos(a;) 003(2) 
— 003(2/) 003(2) 
003(2/) cos(x) 




(17) 



r 



where A'' is a normalization factor such that (a*^*^^)^ 
In the continuum limit this reduces to 



IV. STABILITY OF THE MEAN-FIELD 
SOLUTION: THE ROLE OF THE PSG 



v(0) 



J- J- J-' 

^ V3 V3 V3 







-1 

V3 

-1 -1 

1 1 



1 -1 

7f 



^ J_ ^ 

V3 V3 V3 



(18) 



a^"' acts as a spectrum inverting operator on H, inter- 



changing hole states at energy 
at energy E{k). 



-Eik) with particle states 



We point out that many of the interesting features of 
the low-energy spectrum of the monopolc flux state can 
be generalized to a class of lattices whose geometry is 
related to certain representations of Lie groups^^. In- 
deed, the four sites in the tetrahedral unit cell can be 
viewed as the four weights in the fundamental represen- 
tation of SU(4); hopping on the pyrochlore is then analo- 
gous to acting with the appropriate raising and lowering 
operators. This perspective gives an explicit connection 
between the hopping Hamiltonian (jl3p and the ladder 
operators in the fundamental representation of SU(4). 
Analogous hopping problems can be studied for various 
other Lie group representations, as outlined in detail in 
Ref. m. 

To summarize, the monopolc flux state is a spin liq- 
uid which preserves all symmetries of the full Hamilto- 
nian except P and T. At mean field level it has gaplcss 
spinous along a 1-dimcnsional Fermi surface of 4 lines 
which intersect at the origin. Though strictly at = 00 
it has higher energy than the dimerized state, Gutzwiller 
projection suggests that for = 2 this is no longer the 
case, and the monopole flux state is the lowest energy 
simple mean field ansatz after projection. We now turn 
our attention to what can be said about the stability of 
this rather unusual mean field state. 



Next, we would like to address the question of whether 
the mean field solutions described above maintain their 
basic properties at finite N and whether this holds all the 
way to = 2. This is a difficult problem, whose com- 
plete solution is not available even for the longer studied 
cases of the algebraic spin liquids in d = 2 ^ However, 
following that work the general idea would be to try and 
understand if the state is truly stable at large enough 
A^ while leaving the question of stability at small A^ to 
detailed numerical investigation. 

There are several questions here. First, is the mean 
field solution locally stable? Second, is it the global min- 
imum? Third, assuming the answer thus far is in the 
affirmative, is the expansion about the mean field solu- 
tion well behaved? Ideally, this would mean convergent, 
but it would be sufficient to know that it does not destroy 
the qualitative features of the gapless spinon dispersion 
at mean field. For example, in the case of the algebraic 
spin liquids in d = 2 the spinous interact and acquire 
anomalous dimensions away from N = 00 but they re- 
main gapless in the vicinity of a discrete set of points^^. 
Finally, what is the spectrum of collective (gauge) exci- 
tations that arise in this expansion? 

Based on the experience with spin liquids in d = 2, an- 
swering the first two questions in the affirmative is likely 
to require the addition of more terms to the Hamiltonian 
although it may be possible to choose them so that they 
become trivial at A^ = 321. We have not investigated this 
in detail but there does not appear to be an obstacle to 
doing this. 

The third and fourth questions require detailed consid- 
eration of the symmetry properties and the detailed dy- 
namics of the expansion which is that of a lattice gauge 



See Refs. |29||37| and references therein. 
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theory with matter and gauge fields in some fashion. In 
this work we will carry out the first part of this program 
which goes under the study of the "Projective Symme- 
try Group" (PSG) discussed in detail by Wen2^. In this 
section we review the concept of the PSG and its implica- 
tions for perturbative expansions. We also show that at 
N — oo, or in mean field theory, the PSG already helps 
us understand the stability of particular mean field so- 
lutions; to our knowledge this particular aspect has not 
appeared in the literature before. 

Turning first to the PSG, observe that though the orig- 
inal Hamiltonian formulated in terms of spin operators is 
invariant under the full space group of the pyrochlore lat- 
tice, the actual mean field Hamiltonian of the monopole 
flux state is not: many of the symmetry transforma- 
tions map the mean field Hamiltonian into different but 
gauge equivalent Hamiltonians. Thus, when working in 
the gauge theory formulation of the problem, the actual 
symmetry transformations of the mean field Hamiltonian 
have the form: 

c,^gs{s{c,)) (19) 

where s is an element of the space group, and g is a gauge 
transformation. As the full Hamiltonian is gauge invari- 
ant, ()19|1 is simply an alternative formulation of the lat- 
tice symmetries. Hence as emphasized by Ref. [2^, these 
projective symmetry operators are exactly analogous to 
lattice symmetries in the original spin problem. Indeed, 
the correct choice of gauge transformation ensures that 
both Hmf a-nd H are invariant under the PSG, so that 
the family 

Hx = Hmf + KH - Hmf) (20) 

is also invariant and perturbative corrections in H — Hmf 
cannot break the PSG symmetry. 

Before discussing the implications of PSG symmetry 
for the monopole flux state, we would like to briefly un- 
derline how the PSG constrains the mean field theory at 
infinite which is a much simpler but still instructive 
exercise. 

Ignoring the dimerization instability, the monopole 
flux state is a mean field minimum for nearest-neighbor 
couplings. The PSG is the symmetry group of the corre- 
sponding mean field Hamiltonian. We may now ask what 
happens to the PSG if further neighbor couplings are in- 
cluded in the Hamiltonian, in particular do they lead to 
terms in the new mean field Hamiltonians that modify 
the PSG found earlier? 

At A = 00 this is a problem of minimizing the expec- 
tation value of the sum of the quadratic Hamiltonian in 
Eq.(|4]) and the new generic terms 

SH^-J2J2(^tM^,+h.c.) + ^Yl IX^^-I' (21) 

a {ij) 'J <ij> 

wherein the primed sum runs over non-nearest neighbor 
bonds and the J^, are much smaller than the nearest 



neighbor J. We will now show that, gcnerically, the result 
of the new minimization for the perturbed problem pre- 
serves the PSG for the nearest neighbor problem. While 
we use the language of perturbing about the monopole 
flux state, the argument is general. 

With the addition of the perturbation, the functional 
that we need to minimize over the full set of {Xy} is: 

Emf = {H + 5H)h+sh + 7 E I^^^ I' + |^ E Ix^.f (22) 

Let Xij^ denote the values of the link fields when SH = 0, 
i.e. in the monopole flux state. For small Jy wc expect 
the new minimum to lie not far from the old one, whence 

the link fields will be close to the values Xij''- Conse- 
quently we will compute the expectation value required 
in the above equation in perturbation theory in 6H about 
H. (If such an expansion fails to have any radius of con- 
vergence then we are already parked at a phase transition 
and no stability argument is possible.) 

This expansion, 
EmfUx.,}) =Eo + {0\SH\0) + E + • • • 

ri>0 ° " 

+7Eix.i^ + f Eixdt^3) 

where the numerical indices refer to the ground and ex- 
cited states of the unperturbed Hamiltonian i^({Xij'})j 
has three properties that we need. First, the linear term 
takes the explicit form, 

-Y.Y.^''lc.CJc.)Sx^J+h.c. (24) 

" {ij) 

where Sxij is Xij for the new bonds and the deviation 
from Xij^ for the nearest neighbor bonds. This implies 
that new minimization likes to turn on exactly those 
Xij that transform as the expectation values {cl^Cja}- 
If these are, in fact, what get turned on, then the new 
mean field Hamiltonian will indeed inherit the PSG of 
the starting one. The second property that we need can 
be established by considering a decomposition of SH into 
a piece that commutes with the PSG generators and an- 
other piece that docs not. It is straightforward to see that 
terms from quadratic order and beyond must give rise to 
a potential which is even in powers of the non-PSG con- 
serving piece of SH. Finally, at sufficiently small the 
potential for the Xij must be stable due to the explicit 
factors of l/Jij- Together these properties imply that 
the new minimum must be in the "direction" selected by 
the linear term and hence will exhibit the same PSG as 
before. 
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V. THE PSG OF THE MONOPOLE FLUX 
STATE 

We will now describe the PSG of the monopole flux 
state, and its implications for stability at the mean field 
level. The space group of the pyrochlore lattice is Fd3m, 
which contains 24 symmorphic and 24 non-symmorphic 
elements. For our purposes it is most convenient to divide 
these elements into the 24 proper elements composed of 
rotations and translations, and 24 improper elements in- 
volving a reflection or inversion. The 24 proper elements 
are: 



P" = {1,8C3,3C2,6(72,6C4} 
The improper elements consist of 



(25) 



(26) 



where g denotes a non-symmorphic operation, in which 
rotations or reflections arc accompanied by translation 
along an appropriate fraction of a lattice vector. P° is 
a proper subgroup of Fd3m, while is generated by 
the product of the inversion operator (inversion is taken 
about one of the lattice sites) with the elements of P°. 
The symmetry transformations, along with the full action 
of the PSG, are outlined in Appendix [Bl 

The PSG of the monopole flux state has the following 
general structure, outlined in more detail in Appendix 

m 

• Translations : FCC translations, combined with the 
identity gauge transformation. 

•P*' space group elements : These elements are 
symmetries when combined with appropriate gauge 
transformations, which induce a tt phase shift at some 
the sites in the unit cell. 

•P' space group elements : These elements are 
symmetries when combined with an appropriate gauge 
transformation, as above, and a time reversal transfor- 
mation. 

• Charge conjugation C : The charge conjugation 
operator maps — > c|. 



Restricting Perturbative Corrections Using 
PSG Invar iance 



To deduce what restrictions PSG invariance imposes 
on the spectrum, wc begin with a generic 4x4 quadratic 
Hamiltonian 



(2) _ 



ij 



Jin C," Cj 



(27) 



The bonds Jy connect arbitrary sites in the lattice, but 
respect the lattice symmetries. In what follows, we will 
use the PSG to restrict the possible quadratic terms, and 
show that all terms allowed by symmetry vanish at the 



Fermi surface. Hence the Fermi surface of the monopole 
flux state is unaffected by PSG-preserving perturbations 
to the Hamiltonian. For simplicity we will drop the su- 
perscript (2) in the remainder of this section to simplify 
the notation. 

Though the inversion P and time reversal T are broken 
in the mean field state, the combination PT leaves both 
the full and mean field Hamiltonians invariant. Terms 
invariant under this transformation have the form: 



(J' + iJ")clc,+s + (f - iJ")clc.-s 



(28) 



and the Hamiltonian is real in momentum space. Fur- 
ther, invariance under charge conjugation forces all spa- 
tial bonds to be purely imaginary: under C, 

(J' + iJ")clc.+s + (J' - iJ")cl+sCx 
^ (- J' + iJ")cic^+s + {-J' - iJ")cl^sC, (29) 

so that J' = if C symmetry is unbroken. In momentum 
space, if we write the Hamiltonian as t(;H{k)ip, Equations 
(PS)) and imply that H{k) is real and an odd function 
of k. We may express elements of the matrix i?(k) as a 
superposition 



ifa(,(k) = ^ JR;ab sin [k • (R + Tab 



(30) 



R 



where R is an FCC lattice vector, the indices a, b label 
sites within the unit cell, and Juiat is the coupling be- 
tween sites a and b separated by the lattice vector R, and 
the vector Tab in the unit cell. This is the general form 
for a function periodic in the Brillouin zone. 
Diagonal Terms Let H11...H44 be the diagonal ele- 
ments oi H . To restrict the form of Hn, we consider 
the action of all PSG operations that map site 1 in 
the tetrahedral unit cell onto itself. These are (see Ap- 
pendix [B] for labels and actions of the PSG elements) 
{Ci, C^, 6*23, 6*24, C'34}, which transform Hn in the fol- 
lowing way: 

Hiii^xT ky, kz) > Hi\{kz,kx,ky) 

-4 Hii{ky,kz,kx) 

Hii{-kx, -fc^, -ky) 

Hii{ kz, ky, kx) 

Hii{-ky,-kx,-kz) (31) 

which allows us to express Hn in a form where its sym- 
metries are manifest as: 

Hii{kx,ky,kz) = ^[Hii{kx,ky,kz) + Hii{kz,kx,ky) 

+Hii{ky,kz,kx) - Hii{kx,kz,ky) 
-Hii{kz,ky,kx) - Hii{ky,kx,k^'^Y) 

Similarly we can relate i?22 7^^331 and to Hn by 
considering operations which interchange site 1 with sites 



C-2 
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2, 3, and 4 respectively. These imply: 

H22{kx,ky,kz) = Hii{ky,kx,-kz) 
H33{kx,ky,k^) = Hii{k,,-ky,k^) 
H44 {kx, ky, kz) = Hi i{^kx, kz, ky) 

(33) 

While multiple transformations map between each pair 
of diagonal elements, the group structure and invariance 
of Hii under PSG transformations ensures that these 
mappings all yield the same result. 

The reader should note that Eqs. ([5^ and (|33p ensure 
that along the Fermi lines k = ±(1,±1,±1) all allowed 
diagonal terms vanish. 

It is worth digressing to make one more comment on 
the diagonal terms. Using the symmetrized form of Hn 
in Eq. ([32|) above, we can rewrite the term in Eq. ([30]) 
with a fixed R and a = = 1 as 

>/r;11 sin(k • R) = ^JR-ii[sm(kxRx + kyRy + k^Rz) 

+ sm{kzRx + kxRy + kyRz) 
+ sm{kyRx + kzRy + kxRz) 

- s\i\{kxRx + kzRy + kyRz) 

- sui{kzRx + kyRy + kxRz) 

- sm{kyRx + kxRy + kzRz)] (34) 

The form (j34p vanishes if any two coefficients are equal; 
non- vanishing terms occur only for a sum of at least three 



FCC translations. Physically this corresponds to a hop- 
ping between a site and its translate some three lattice 
vectors distant. 

OfF-Diagonal Terms As H is real in momentum space, 
Hab = Hba- To restrict the form of H12, consider the 
action of all PSG elements which either map sites 1 
and 2 to themselves, or interchange them. These are 
{C34, C12, Cz}, which transform H12 according to 

Hi2{kx,ky,kz) ^* -Hi2{-ky,-kx,-kz) 
^ H2i{ky,kx,~kz) 
^ -H2i{-kx,^ky,kz) (35) 

Again, transformations mapping sites 1 and 2 onto other 
sites in the unit cell can be used to deduce the form of 
the remaining off-diagonal elements. Hence 

Hi2{kx, ky, kz 

kx ; ky ) 

Hi4 {ky,kz,kx) 

-H23{-ky,kz,-kx) 

S -H2i{-kz,kx,~ky) 
c 

-H34{-kx,ky,-kz) (36) 
This gives off-diagonal entries: 





Hi2{kx, ky, kz 



Hi2{kx,ky,kz) Hi2{kz, 



Hi2{kz,kx,ky) Hi2{ky, 
Hi2{ky,kz,kx) i?i2(fczr 



Hi2{ky,-kz 
,kx) 



Hi2{ky,kz,kx) 
kx) Hi2{kz,-kx,ky) 
Hi2{kx,-ky,kz) 



kxjky) Hi2{kx, ky,kz) 







(37) 



r 



where again we can make the symmetries manifest by 
writing 

Hi2{kx,ky,kz) = ^[Hi2{kx,ky,kz) + Hi2{ky,kx,kz) (38) 

+Hi2{ky,kx, ~kz) + Hi2{kx, ky, -kz)] . 

Again, it is useful to focus on the contribution to H12 
from bonds with a given R which can now be seen to 
come with the factor: 

cos(fc2i?j,)[sin(fc^i?^ + kyRy + 

+ sm{kyRx + kxRy + ^l±hL)] . (39) 

Eq. ([M)) shows that Hi2{kx,ky,kz) vanishes along the 
lines {k, —k, —k), {—k, k, — fc), (— fc, — fc, k). Of course, this 
can also be seen directly from Eq. ([55)) . 



Now we may consider the fate of the monopole flux 
state's exotic Fermi surface. Since PSG rotations map 
between different Fermi lines, it is sufficient to con- 
sider possible alterations to the spectrum on Fermi line 
[k, k, k). The most general form that H can have about 
the line {k, k, k) is: 

Hi2{k,k,k) Hi2(k,k,k) Hi2{k,k,k) 

^_ Hi2{k,k,k) 

Hi2{k,k,k) 

Hi2{k,k,k) 

(40)' 

which has two zero eigenvalues. Thus terms allowed by 
symmetry add neither a chemical potential nor a gap 
to any part of the Fermi lines, and preserve the charac- 
teristic structure of the monopole flux state, with 2 low 
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energy states about each Fermi line, and 4 low energy being modified as a function of the momentum along the 
states about the origin. line. Indeed, (|l(7)l implies that the general form of vp is: 

Note that nothing prevents the Fermi velocity vp from 



vpik) 



= jR;i2(i?x + Ry + 1/2) cos(2fci?^) cos2k{Ry - R^) + ^^^^{Ra, sin(2fci?^) sin2fc(i?^ - Ry) 



-Ry sm{2kRy) sin 2k{R^ ~ R^) + R^ sm{2kR^) sin 2k{Ry - Rx)) 

I 



(41) 



B. Time reversal and Parity 

One striking feature of Hmf is that it is odd under 
both T and P, reminiscent of the chiral spin state first 
described in Ref. [H. Though T is naively broken, some 
care must be taken to show that the apparent T break- 
ing is physical and that \il'),T\'ip) are gauge inequivalent 
states^. Readers familiar with this subtlety from discus- 
sions of T breaking on the square lattice, should note that 
the pyrochlore lattice is not bipartite and hence naive 
time reversal is no longer equivalent to particle-hole con- 
jugation. But most directly, as explained in Ref. |23| . the 
operator 



tries of H 



MF 



Eijk — Si ■ {Sj X Sk) 



(42) 



where the spins i, j, and k lie in a triangular plaquette, 
is odd under T and P. Hence if {Eijk)\^) ^ 0, the state 
breaks time reversal. 

At mean field level, 

— ^ 

^(-Bl23) = (Xl2X23X3l) - (Xl3X32X2l) (43) 

and states with an imaginary fiux through triangular pla- 
quettes are T-breaking. For the monopole flux state, we 
have confirmed numerically that this T-breaking is ro- 
bust to Gutzwiller projection; the results are shown in 
Table in 



Lattice size 




3x3x3 
5x5x5 


0.039 
0.043 



TABLE XL Expectation values of the T-breaking operator 
(Ba) for triangular faces of the tetrahedra. 



We also note the curiosity that at infinite N , the 
spectrum-preserving nature of T and P allows us to 
construct additional symmetries which are noi, however, 
symmetries of the full H . Particle-hole symmetry at each 
k allows us to construct the following 2 discrete symme- 



f : |V;(x,t)) ^aW|^(x,-i)) 
T: |i^(rr,i)) ^a(°)|^(-x,t)) 



(44) 



where was defined in Eq. pT|) . Both of these commute 
with the non-interacting Hamiltonian: since aQ^Ha^ = 
— i?, we have 

(V'lTtHfiv) = 

= (V'|i?|^) 

{i}}\P^HP\^) = 

= (45) 

The matrix structure of ao is such that T and P are not 
symmetries of the full Hamiltonian, however, and will not 
be robust to perturbative corrections about mean field. 



C. PSG symmetry and perturbation theory in the 
long wavelength limit 

We have established that invariance under the PSG 
transformations and charge conjugation forbid both mass 
and chemical potential terms on the Fermi lines. Here 
we explore how these PSG symmetries are realized as 
symmetries of the linearized low-energy theory away from 
the origin, and hence see in that setting why they are 
protected perturbatively. 

Consider the linearized theory about the Fermi line 
Zi = (fc, k, k). A general Hamiltonian in the 2x2 space 
of low-energy states can be expressed as: 

i?(k) = Vl(k)/i(k)V'i(k) (46) 
/i(k) = ^(fc) + TO(fc)fT3 + £(fc, 'i;)(cos(0)cri -f- sin(6')(T2) 

where k is the component of the momentum k along the 
line, and (u, 6) are the magnitude and angle respectively 
of the momentum perpendicular to the line. Here 



V'ii(k) 
V'i2(k) 



(0,W^1,LJ) 
(0,C^,1,C^2) 



.(k) 
.(k) 



(47) 
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with Lo = e^'^*/^. The states (|47|) are eigenstates of the 
rotation operator Ci which rotates about the (1,1,1) di- 
rection, with Ciipij = Lu^ipij. 
Under charge conjugation. 



Vii(k)^Vl2(-k) 
Vi2(k)^Vli(-k) 



(48) 



The corresponding symmetry operator in the continuum 
theory is 



C:^/'i(«,0)^c7i[Vl(«,7r-0)]^ 



(49) 



with the Fermi surface points at k and — /c interchanged. 
This imphes m{—k) = m{k), and n{—k) = —i.i{k). 

Further, an analysis of the PSG transformations re- 
veals that the glide rotations (C^) map clockwise rotat- 
ing states to counter-clockwise states while reversing the 
direction of the corresponding Fermi line: ipn^kh) — > 
^i2{—kli). This transformation leaves the mean field 
Hamiltonian invariant. In the 2x2 basis, this is because 



T: |V(i',e)) ^al|V(^^7^-^^)) 



(50) 



is a symmetry of the mean field Hamiltonian. Note that 
the momentum transformation can be realized in 3 di- 
mensions by a TT rotation about the line x = y, and hence 
should also send k — fc, though there is no way to de- 
duce this from the form of the mean field eigenstates. 
The symmetry transformation ([50]) reverses the sign of 
the mass term, but not of the chemical potential, im- 
plying that m{—k) = —m{k) and /i(fc) = /x(— fc). Hence 
we conclude that in the continuum theory about a given 
Fermi line, the symmetries C and T prevent a gap or 
chemical potential from arising. 

One might ask why we have not considered mass gaps 
of the form mai or m<T2 ; both of these choices turn out to 
violate either or (|50p . Indeed, both choices explicitly 
break the rotational symmetry of the spectrum about the 
Fermi line. 



VI. CONCLUDING REMARKS 

We have discussed an interesting mean field (large N) 
solution to the Heisenberg model on the pyrochlore lat- 
tice. This is a P and T breaking state in which all tri- 
angular plaquettes have an outward flux of 7r/2. After 
Gutzwiller projection, this state has lower energy than all 
other mean field states considered, including the simplest 
dimerizcd state. Its low-energy physics is rather striking, 
with a spinon Fermi surface of lines of nodes preserving 
the discrete rotational symmetries of the lattice. The 
symmetries of the Hamiltonian suggest that this Fermi 
surface is perturbatively stable and thus should charac- 
terize a stable spin liquid phase, at least at sufficiently 
large N . 

However, our analysis of stability is thus far based 
on only on symmetries and does not rule out dynami- 
cal instabilities. The study of such instabilities requires 



adding back in the gauge fluctuations that are suppressed 
at = oo and studying the coupled system consisting 
of spinous and gauge fields. In the well studied case of 
two dimensional algebraic spin liquids it took a while to 
understand that this coupled system could, in fact, sup- 
port a gapless phase at sufficiently large N despite the 
compactness of the gauge fields. In the present prob- 
lem there is also the specific feature that at A = 2, as 
the background flux per plaquette is t^(l), such an anal- 
ysis should incorporate fluctuations of an SU{2) gauge 
field2^.We defer addressing this set of questions to future 
workSl. 

Finally, we note that our initial motivation in this 
study was to see if we could construct a fully symmetric 
spin liquid on the pyrochlore lattice for S' = 1/2 in con- 
tradiction with previous studies using other techniques. 
We have not succeeded in that goal and, as the technique 
in this paper has produced a pattern of symmetry break- 
ing distinct from the ones considered previously, the fate 
of the S = 1/2 nearest- neighbor Heisenberg antiferro- 
magnet on the pyrochlore lattice remains undeciphered. 
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APPENDIX A: MONTE CARLO AND 
PROJECTED WAVEFUNCTIONS 

Gutzwiller projection can be carried out exactly for fi- 
nite systems using the projector in the Slater determinant 
basis, using the method of HI- At half filling, the Slater 
determinant is represented as a product of a spin up and 
a spin down determinant of equal size {N/2) where N 
is the number of sites in the lattice. Each site is rep- 
resented exactly once in either the spin up or the spin 
down matrix, yielding a wavefunction which obeys the 
single occupation constraint on all lattice sites and is a 
total spin singlet. The problem then reduces to the evalu- 
ation of expectation values of operators for wavefunctions 
of finite systems with definite spin distributions on the 
lattice: 



(0) = 



{i'\o\i') 
ii'H') 



(Al) 



where \ip) = '^ai'^Wl'^)' ^^'^ 1*^) ^ specific distribu- 
tion of spins on the lattice. 



H=n4.Tn4,aio)- 



(A2) 



The expectation value of the operator O is evaluated by 
summing over all spin configurations on the lattice. To 
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evaluate this sum we follow the approach of |3J, which we 
will review here. The expectation value is given by: 



(O) 



where 



a 13 



{a\o\p){m\\{c^m 



OL 

(A3) 



/(«) 
p{a) 



a p 



(A4a) 
(A4b) 



Here p{a) > and J2aP('^) ~ 1' which makes p{a) a 
probability distribution. The expectation value can be 
rewritten to resemble a weighted sample using p{a) eval- 
uated using a Monte Carlo sampling. The evaluation 
is executed using a random walk in configuration space 
with the weight p{a). The transition probability Taa' of 
the Monte Carlo step is 



p{a')/p{c 



p{a') > p{a), 
p{a') < p{a). 



(A5) 



The configuration a' is generated by exchanging a ran- 
domly selected pair of oppositely oriented spins. We also 
calculate various operators in the mean field wavefunc- 
tion to test the accuracy of the algorithm. In this case 
the one particle per site constraint is not imposed; the 
configuration a' is generated by moving an up or down 
electron at random to another empty up or down site. 

Pyrochlore has an FCC lattice with a four point basis. 
We use the rhombohedral unit cell for the Monte Carlo 
evaluation, using boundary conditions periodic along the 
fee directions. The spin correlations turn out to be quite 
insensitive to boundary conditions for the lattice sizes 
that we have considered. 

To compute the Slater determinant, the mean field 
Hamiltonian (jTU]) is diagonalized in the band eigenstates 

H = ^ cl%Hab{k)cbka = ^ E^{k)al'la„ka, (A6) 



„6,fc 



where cj^^ is the Fourier transform 



4" _ ik-{R,+raj2) Ja 



/2- 



(A7) 



Here Tai/2 refer to the points in the four site basis of 
the tetrahedral unit cell. If we assume that the bands 
J/ = 1, 2 are filled, the operator cj^^ can be expressed as 



(A8) 



where the matrix S\ ^(k) diagonalizes the Hamiltonian. 
The mean field ground state is just the Fermi sea filled 



to the appropriate Fermi level which, in this case, is the 
1st Brillouin zone boundary. 

m,nean= J] «U4l|0) (A9) 
k<kf ,a 

We rewrite the mean field Fermi sea |$)meara m first quan- 
tized form; 



U,D 



(AlO) 



where 



^ ■■■'■^./.. tI0)' (All) 



(A12) 



The basic Slater determinant wavefunction for each spin 
orientation is of the form 



^{Ri,ai,kj,v)a = Det[(/)(fii,ai : kj,v) 



(A13) 



where (j){Ri^ai '■ kj,v) is the single particle wavefunc- 
tion of the electron at Ri,ai = Ri + I'ai/'^- The wave 
number kj refers to a point in the conjugate lattice in 
the rhombohedral Brillouin zone, and v is the band in- 
dex. In terms of the matrix Sa.v{k), the single particle 
wavefunction is 



4'(.Ri,ai ■ kj , I') — Sa.u{kj] 



./2) 



(A14) 



In the above, U = {Ri,ai, -- ^ RNL/2,aN^/2} is the 
set of lattice sites occupied by up spin electrons, and 
D = „ , • • • , R'r.T /n I is the set of lattice sites 

occupied by down spin electrons. 

Gutzwiller projection is imposed by ensuring that the 
two sets U and D have no elements in common. The 
Monte Carlo update is performed by exchanging rows 
selected at random from the $[/| and $d| matrices. 

The calculation of the transition probability Taa> in- 
volves the calculation of determinants of matrices of size 
Nl/2 X A^l/2, an 0(7V|) operation. The algorithm of 
Ceperley, Chester and Kalos^ reduces this to an 0(iV|) 
operation for the special case of updates involving one 
row or column. The matrix = ^jji a-nd the trans- 
pose of its inverse (similarly for $ di ) arc stored at 
the beginning of the Monte Carlo evolution. If the up- 
date changes the ath row Aj,M^ -> M'^ , the 
transition matrix is T^a' = DetM'^/DetM^: 



/u 



DetAf 
DetM^ 



(A15) 



This is due to the fact that Dct(M'^A'/ ) is the matrix 
of cofactors. If the move is accepted, the inverse matrix 
can be updated using 0(iV|) operations: 



M 



I = a, 

Ml - Ef=/' A,Ml/r ^ + a. 



(A16) 
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The evaluation of operator expectations values has to 
be handled with care as we are dealing with fermions. 
The relative sign of determinants must be tracked in a 
consistent way. Thus we write all spin configurations in 
the order: 

la) = ct, t'^'ctj ticL, |---ct,, I 10). 

(A17) 

The wavcfunction is given by ($|q!) — <E>c/| $dx (o;) , the 
Slater determinant eigenfunction with the given spin dis- 
tribution a. We are interested in the operators ^ S'|j. 

^R, „ . and ^.-Sn^^ +S-^^ -Sl, ). To keep 

track of the proper sign specification we express the spin 
operators in terms of fermionic operators with the same 
order as the spin configuration: 

(A18) 

The amplitude '^r , I*^) the determinant wave- 

fimction <i>c/|$D| (a') with the rows of <i>[/| and ^di 
changed as described above. S'|j. • \a) is easier 

to evaluate as it is diagonal. 

As a result of the four site basis there are only N^/A 
lattice points in the Brillouin zone. Therefore, the ac- 
curacy of the Monte Carlo evaluation is limited by the 
number of points in k space. We have used a lattices of 
size 5 x 5 x 5 (500 sites). 

A basic Monte Carlo 'move' consists of 2Nl updates 
followed by a sampling. The first 10000 moves were used 
for thermal relaxation and were discarded. 50000 sam- 
ples were used for the evaluation of expectation values. 
These were divided into 10 sets and the average in each 
set was used to estimate the statistical fluctuations and 
error bars of the Monte Carlo evaluation. To check the 
accuracy of the algorithm and the effect of finite lattice 
size we evaluated spin correlation functions of the mean 
field states and compared with results from the numer- 
ical evaluation of Green's functions. The Monte Carlo 
results are quite close to the expected values (see Table 
imp . The site indices of the spins refer to Fig. [H 



APPENDIX B: THE PSG OF THE MONOPOLE 
FLUX STATE 

1. Symmetries of the Pyrochlore Lattice 

The space group FdSm of the pyrochlore lattice con- 
sists of the 24 element tetrahedral point group 43m, and 
a further 24 non-symmorphic elements. We will briefly 
describe the actions of these symmetry operations here. 



TABLE III: Spin Correlations from Variational Monte Carlo. 
To check the accuracy of the algorithm, correlators are also 
calculated at mean-field, using both Monte Carlo (MC) and 
Green's function (G) approaches. 




FIG. 4: Sites for the spin correlation functions calculated in 
Table Uni 



All vectors are expressed in the basis of the standard 
cubic FCC unit cell. 

The actions of the tetrahedral point group fix the posi- 
tion of one tetrahedron's center (at e.g. (a/8, a/8, a/8)). 
Its elements are: 

1. The identity 

2. 8C3 : there are four 3-fold axes, one passing 
through each vertex and the center of the oppo- 
site face. Rotations about this axis permute the 
3 vertices not on the axis. These we label C1..C4, 
Cf..C^, where Ci,Cf fix site i of the tetrahedral 
unit cell. 

3. 3C2: There are 3 2-fold axes, parallel to the x,y, 
and z axes. Each axis bisects a pair of edges on the 
tetrahedron; the ensuing rotation exchanges pairs 
of vertices. These 

4. Gfjf/: A plane of reflection passes through each edge, 



Trial Wavefn. 


S!5| 


S'l 5*3 


SI Si 


SI SI 


mean field 










Flux (G) 


-0.01388 


-0.00188 


-0.00097 


-0.00097 


Flux (MC) 


-0.01386 


-0.00192 


-0.00103 


-0.00084 




±0.00003 


±0.00004 


±0.00005 


±0.00006 


Monopole (G) 


-0.01745 


0.00000 


0.00000 


-0.00087 


Monopole (MC) 


-0.01713 


-0.00002 


-0.00008 


-0.00085 




±0.00004 


±0.00001 


±0.00002 


±0.00006 


Projected 










Flux 


-0.04169 


-0.00137 


0.0029 


-0.00270 




±0.00004 


±0.00005 


±0.0001 


±0.00008 


Monopole 


-0.0497 


0.00631 


0.00528 


-0.00499 




±0.0001 


±0.00002 


±0.00004 


±0.00005 
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and out the center of the opposite face. These 
planes lie on the diagonals with respect to the FCC 
cubic unit cell, and hence are called diagonal reflec- 
tions. 

5. 654: An improper rotation of degree 4 about the 
axis bisecting 2 edges (parallel to the x, y, or z axis) 
is also a symmetry. The tetrahedron is rotated by 
7r/2 about e.g. (1/8,1/8, z) and reflected through 
the plane z = 1/8. This operation squared pro- 
duces one of the 2-fold rotations, so each axis con- 
tributes 2 group elements. 

The remaining non-symniorphic elements, which we 
distinguish from their symmorphic counterparts using 
the notation S, are: 

1. 6C4: There are three 4-fold screw axes: 
(3/8, 1/8, z), (3/8, y, 1/8, and (x, 3/8, 1/8). The 
symmetry rotates the lattice by 7r/4 about such an 
axis, and translates by 1/4 of the side length of 
the FCC cubic unit cell along the axis. Each axis 
accounts for 2 elements of the quotient group, as 
C4 = tC2, with t an FCC translation and C2 one 
of the 2-fold rotations of the point group. These we 
label Cx , C,j. , Cy , Cy , Cz , C^ . 

2. 6C2' Along each of the 6 edges of the tetrahedron 
(the FCC basis vectors) there is a 2-fold screw axis. 
The lattice is translated along the edge of a tetra- 
hedron, then rotated by n about this edge. These 
we label Cy , where Cij has a screw axis along the 
line joining sites i and j. 

3. 3<7h- The x,y and z planes of the cubic unit cell 
each contain a horizontal glide plane. The lattice 
is translated along an FCC vector in the plane, 
e.g. by (1/4,1/4,0), and then reflected through 
the plane - in our example, through z ~ 0. 

4. i The lattice is inverted about the origin. 

5. 856: The products of the 8 C3 rotations with the 
inversion give 8 improper 3-fold rotations. (These 
are not in the point group because they map a sin- 
gle tetrahedron onto its neighbor.) 

For our purposes these 48 elements divide into 24 C el- 
ements involving pure rotations and translations, and 24 
S elements involving improper rotations, reflections, or 
inversions. The S elements are not symmetries of the 
monopole flux state, as they map monopoles to anti- 
monopoles; to construct the appropriate symmetry ele- 
ments they must be combined with a time reversal trans- 
formation. Since all such elements can be expressed as a 
product of a C element with the inversion, this is simply 
a consequence of the fact that while P and T are sepa- 
rately broken in the monopole flux state, the combination 
FT is still a symmetry. 



2. The PSG 



Here we list the PSG transformation rules for the sym- 
metry operations described above. Throughout, we use 
the gauge illustrated in Fig. [U in which all bonds are 
either ingoing or outgoing from site '1'; starting from a 
different gauge will permute the gauge transformations 
listed here (note that this has no elfect on which bonds 
are allowed or disallowed by the PSG, however). Note 
that these tables only show the mapping between the 
site labels 1...4; it is important to bear in mind the ef- 
fect of the translations in the case of the non-symmorphic 
elements, which reverse the directions between sites by 
interchanging up and down triangles. To this end we 
also include a table of momentum transformations under 
these operations. 

Since S elements are products of C elements and inver- 
sion, it is sufficient to consider the 24 rotation operations, 
together with the operator iT. Note that group multi- 
plication in the PSG is valid modulo the global gauge 
transformation Ci — > — q, which clearly does not alter 
the Hamiltonian. Thus only the relative phases of the 4 
sites in the tctrahedral unit cell are relevant to the PSG 
transformation. 



TABLE IV: Action of PSG Point Group Rotations on spinon 
operators. 



1 


Cl 


Cl 


C2 


Cl 


Cz 


Ci 


C4 


Cl 


c. 


Cy 


Cz 


Cl 


Cl 


Cl 


Ci 


-C3 


C2 




C3 


— C2 


— C4 


-C3 


C2 


C2 


C3 


C4 


— C2 


— C2 


C4 


Cl 


— Cl 


C3 


-C3 


C4 


— Cl 


C3 


C4 


C2 


-Cl 


C4 


-C3 


-C3 


C2 


Cl 


C2 


Cl 


C4 


C4 


C2 


C3 


C3 


Cl 


— Cl 


C2 


— C4 


— C4 


Cl 


— C2 


-C3 






ky 




ky 


kz 


— ky 


-kz 


— ky 


kx 


kx 


kx 


ky 


kx 


kz 


kx 


-kz 


~~kx 


-kz 


kx 


kz 


— ky 


ky 


— ky 




ky 


kx 


~ky 


kx 


— ky 


kx 


ky 


kx 


kz 


-kz 


kz 



TABLE V: PSG action of screw rotations 





C12 


Cl3 


Cl4 


C23 


C24 


C34 


Cx 


Cy 


Cz 


Cx 


Cy 


Ci 


Cl 


C2 


C3 


C4 


-Cl 


-Cl 


-Cl 


C2 


C4 


C3 


C3 


C2 


C4 


C2 


Cl 


-C2 


C2 


C3 


C4 


C2 


-C4 


Cl 


C4 


Cl 


C3 


-C3 


C3 


C3 


Cl 


-C3 


C2 


C3 


C4 


Cl 


C2 


— C2 


C4 


-C4 


Cl 


C4 


-C4 


C4 


Cl 


C4 


C2 


C3 


C3 


-C3 


Cl 


-C2 


Cl 


C2 


kx 


ky 


kz 


kx 


kx 


kz 


— ky 


kx 


kx 


kz 


-kz 


— ky 


ky 


ky 


kx 


— ky 


kz 


-kz 


— ky 


kx 


-kz 


kz 


ky 


ky 


kx 


kx 


kz 


-kz 


kx 


ky 


— ky 


kx 


-kz 


ky 


— ky 


kx 


kx 


kz 


kz 
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TABLE VI: Effect of Point Group rotations on low-energy 
eigenstates. 



3. Action of the PSG on Low-Energy Eigenstates 







W2 


Co 


K^4 


n 




n 








— LU^4l 


— LU^'tp21 


—"041 




LJ'021 


1pl2 


UJ^ll'12 


ljJtp32 


— a;^'042 


— LJ'?/?22 


—'042 


Up1p'i2 


a-'^^22 


1p21 


V'41 




^11 






-LJ^t/)41 




1p22 


Tp42 


io'^1p22 


tpl2 


UJ^i'32 


-'032 


-U}i>41 


-1^-012 


1p31 


-Uji>21 


l/'41 




-l/'ll 


i>21 




(^^i>41 


1p32 


—UJ^ip22 


■Ip42 


u'^'4'32 


-'012 


tp22 


-a;V'i2 


1^-042 


tpil 








L01p41 


^11 


ijJ1p2l 


-1^-031 


1p42 


-UJTp32 


UJ^tpl2 


—L01p22 


Up 1p 42 


V'12 


Lo'^1p22 


— a;^032 



Recall that the eigenstates of the low-energy excita- 
tions about the Fermi surface can be expressed in terms 
of eigenstates of the point group rotation operators 8C3 

(c.f. m)- 



TABLE VIL Effect of 2-fold glide rotations on low-energy 
eigenstates. 





C12 


Cv3 


Cl4 


C23 


C24 


C34 


■011 


tJ'022 


LJ^032 


— a;0;42 


-UJ^1pl2 


-'012 


UJ1pl2 


■012 


a;^V2i 


a;03i 


-t^^^41 


-ut/jii 




CJ^011 


'021 


'012 


a;022 


^22 


LJ032 


-'042 


— LJ^022 


'022 


'011 


LJ^021 


^21 




—^41 


— a;022 


'031 


— aj'032 


012 


a;032 


-022 


0)32 


— LJ^042 


'032 




011 


CJ^V31 


-021 


0'31 




V'41 


—UJ^ll'42 


LO^l>42 


-a;^Vi2 


t^'042 


022 


0J1p32 


-042 


-0J1p41 


tJ^041 




0)^0)41 


V'21 


tJ^031 



Tables IB 31 and IB 31 list the action of the rotation el- 
ements of the PSG on these low-energy states. Since 
(|47)) is PT invariant, it is sufficient to consider the action 
of the rotation subgroup; the other PSG elements are 
combinations of rotations with the inversion, and cannot 
procure new information about the low-energy behavior. 
Note that proper rotations always map clockwise rotat- 
ing states to clockwise states, and counter-clockwise to 
counter-clockwise. The improper rotations, conversely, 
map counter-clockwise rotating states into clockwise ro- 
tating states, and vice versa. 
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